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ABSTRACT 

Recent observations have revealed a superbubble associated with the young stellar cluster OC1 
352 near the W4 H II region: a void in H I emission (Normandeau, Taylor, & Dewdney), and a 
bright shell in Ha emission (Dennison, Topasna, & Simonetti). We investigate the hypothesis that 
the bubble is blown by stellar winds from the O-type stars in the association. The Kompaneets 
approximation is adapted to model a wind-blown bubble in a stratified interstellar medium. 
We describe some general principles necessary for understanding the dynamics of an expanding 
bubble and the associated ionization structure in a stratified atmosphere. The Kompaneets model 
can be used to determine the mean scale height of the ambient medium as well as the age of the 
bubble. The ionization structure also places constraints on the ambient density near the cluster. 
We also estimate the surface brightness of the shell and the fraction of ionizing photons which 
escape the bubble. The prescription we use can be applied to any observed bubble that is blown 
by the effectively continuous energy output of stellar winds or multiple supernovae. Application 
to the W4 superbubble shows that the mean scale height of the ambient gas around the cluster 
is remarkably small, 25 pc for a cluster distance of 2.35 kpc. The age of the bubble is estimated 
to be about 2.5 Myr, consistent with the notion that the bubble is blown by stellar winds from 
a very young cluster in which no supernovae have yet occurred. 



Subject headings: H II regions - ISM: bubbles - ISM: individual (W4 supershell) - shock waves 
supernova remnants 



1. Introduction 

Normandeau, Taylor, & Dewdney (1996, 1997; hereafter NTD) have observed a large cavity in H I above 
the W4 H II region in the Pilot project of the arcminute resolution Canadian Galactic Plane Survey (CGPS). 
The conical cavity, found in the velocity range i>lsr — —31.9 to —46.7 km s , opens upward above (and 
perpendicular to) the Galactic plane, and is bounded below by the W4 H II region, which is part of the 
larger W3/W4/W5 H II region complex located in the Perseus arm of the outer Galaxy. Panoramic images 
of CO and thermal dust (IRAS) emission from this region, also at arcminute resolution, have been presented 
by Heyer & Terebey (1998). An association of nine O-type stars, part of the stellar cluster OC1 352 (also IC 
1805), centered at (/, b) = (134.7°, 0.9°), is situated near the base of the cavity, suggesting that the cavity has 
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been created by outflows from these stars. Spectroscopic parallax (Massey, Johnson, & DeGioia-Eastwood 
1995) places OC1 352 at a distance 2.35 kpc; at this distance, the H I cavity has a width ~ 110 pc at a latitude 
b ~ 3.5°, or ~ 110 pc above the cluster. Radio continuum data presented by NTD show the presence of a 
shell of ionized gas along the lower edge of the H I cavity. Dennison, Topasna, & Simonetti (1997; hereafter 
DTS) have observed the same region in Ha with a field of view extending to higher Galactic latitude. Their 
image shows that the H I cavity of NTD is located near the base of a highly elongated shell of Ha emission 
which appears to close far above the Galactic plane. This shell is interpreted here to be the dense ionized 
wall of swept-up gas surrounding a bubble of hot rarefied gas created by the outflow from the massive stars. 

This structure is one example of a vast collection of bubbles and associated shells that are thought to 
populate the interstellar medium (ISM) of galaxies. In our Galaxy, H I observations (Heilcs 1979, 1984) reveal 
numerous shells with dimensions of several hundred parsecs. Other spiral, irregular, and dwarf galaxies also 
show evidence for shells and holes in H I (e.g., Brinks & Bajaja 1986; Deul & den Hartog 1990; Puche et al. 
1992). The holes in H I may be sites where the bubbles have blown out of the stratified ISM and into the 
halo of the host galaxies. Norman & Ikeuchi (1989) suggest that these "chimneys" are a common feature in 
our Galaxy. 

The well-resolved structure of the W4 superbubble and the well known properties of the stellar cluster 
OC1 352 that drives the bubble afford us a chance to apply theoretical models to gain insight into the massive 
star-ISM interaction in a star-forming region of our Galaxy. We show that the highly elongated nature of the 
superbubble implies the presence of significant vertical stratification in the ambient ISM; the scale height is 
remarkably small. 

Theoretical models of bubble expansion in a stratified medium come in many varieties. Kompaneets 
(1960) found a semianalytic solution for the expansion of a blast wave in an exponential atmosphere. The 
Kompaneets approximation assumes the following: 1) uniform pressure within the bubble, 2) bubble expan- 
sion in a direction normal to the local surface, and 3) expansion speed implied by a strong shock (i.e., the 
internal pressure dominates the external pressure). Kompaneets found an analytic expression for the shape 
of the bubble during its expansion. The time-evolution of the bubble can be found by numerical integration. 
The Kompaneets solution has been used and adapted for application to various astrophysical phenomena, 
including rclativistic blast waves (Shapiro 1979), active galaxy winds (Schiano 1985), and impacts within the 
deep gaseous envelopes of giant planets (Korycansky 1992). More sophisticated models of bubble expansion 
include the thin-shell approximation (MacLow & McCray 1988; Bisnovatyi-Kogan, Blinnikov, & Silich 1989), 
which employs the first two assumptions of the Kompaneets approximation, but determines the expansion 
speed through direct numerical integration of the momentum equation for various segments of the thin shell 
of swept-up gas. This numerical approach accounts for the inertia of the swept-up shell, and external pressure 
and gravity can also be included. Finally, full numerical integration of the hydrodynamic equations (e.g., 
Tomisaka & Ikeuchi 1986; MacLow, McCray, & Norman 1989; Tenorio-Tagle, Rozyczka, & Bodenheimer 
1990) and magnetohydrodynamic equations (Tomisaka 1992, 1998) yield the most complete solutions to 
date. For a review of models of bubble expansion in the interstellar medium, see Bisnovatyi-Kogan & Silich 
(1995). 

In this paper, we choose to apply the Kompaneets model because, in addition to its simplicity, the 
semianalytic solution allows a more direct insight into the physics of bubble expansion. An observed bubble 
can be matched to theory using analytic expressions, yielding straightforward estimates of the atmospheric 
structure and bubble age. The original Kompaneets model is adapted so that the energy input is continuous 
rather than occurring at an initial instant, as is expected with stellar winds or multiple sequential supernovae. 
We discuss the conditions under which the Kompaneets solution is valid, and show how some of its properties 



- 3 - 



can be used in an analysis of observed bubbles. Application of the model to the W4 superbubble gives us 
a more accurate age estimate than is possible with models of spherical expansion, as well as a dynamical 
estimate of the gas scale height in the W4 star-forming region. Given the structure of the bubble and ambient 
atmosphere, we also illustrate the evolution of the ionized gas structure created by the O-type stars. The 
current structure of the ionized shell places an additional constraint on the mean density of the ambient 
atmosphere and allows predictions of how much ionizing flux escapes through the bubble. 

The outline of this paper is as follows. In § 2 we apply the Kompanccts model to the structure and 
evolution of the wind-blown bubble. The complementary role of the ionizing UV flux from the massive stars 
is discussed in § 3. We interpret our results in § 4, and a summary of our main conclusions is given in § 
5. Readers who wish to review a more complete exposition of the Kompanccts solution are referred to the 
Appendix. 



The early evolution of a superbubble in a stratified medium can be described by the model for spherical 
expansion of a stellar-wind bubble in a homogeneous medium, as developed by Castor, McCray, & Weaver 
(1975) and Weaver et al. (1977). Their model can be applied to bubbles powered by the combined winds 
of several O and B stars if they are closely associated, and can even apply to bubbles powered by multiple 
sequential supernovae from such an association, since the energy output is effectively continuous (McCray & 
Kafatos 1987; MacLow & McCray 1988). We briefly describe the model of Castor et al. (1975) and Weaver 
et al. (1977) before proceeding to discuss the effects of vertical stratification. 

The earliest phase is one of free expansion of the wind, which drives a shock front into the ambient 
ISM. A shell of swept-up ISM is formed behind the shock front, and once its inertia slows down the shell 
appreciably, an inner shock front is also formed. This inner shock slows down the hypersonic free wind 
before it encounters the more slowly moving shell, creating a region of hot, shocked wind gas. Hence, most 
of the bubble evolution is characterized by a four-zone structure: (a) an innermost freely expanding stellar 
wind surrounded by the inner shock surface, (b) a hot, almost isobaric region consisting of shocked wind 
gas and a small fraction of swept-up ISM gas, (c) a shell of swept-up ISM gas separated from region (b) by 
a contact discontinuity, and (d) the undisturbed ISM, separated from region (c) by the outer shock front. 
Radiative cooling causes region (c) to collapse into a thin, dense shell within a few thousand years, while 
cooling remains relatively unimportant in the rarefied region (b). The inner shock front radius R\ is found 
to be much smaller than the outer shock front radius i?2 • Hence, the expansion of the bubble can be modeled 
effectively as that of a thin-shell driven into the ambient medium with uniform density po by an isobaric 
bubble of radius R s = i?2- The internal energy of the bubble i? t h (and hence the mean pressure P) is 
determined by the energy input from the luminosity Lq of the wind minus the work done in expanding the 
bubble. Under these conditions, Castor ct al. (1975) find the analytic solution 



2. Dynamical Evolution of a Superbubble 



2.1. Early Evolution and Blowout Condition 
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where ft = AttR^/3 is the volume of the bubble, and a ratio of specific heats 7 = 5/3 has been used for 
the equation of state in the bubble interior. The expansion will dissipate and the shell will begin to merge 
with the interstellar medium when the expansion speed u = |i? s A (oc i~ 2 / 5 ) decreases to the sound speed 
of the external medium c SjC — (P c / po) 1 ^ 2 , where P e is the external pressure (essentially the same time that 
P = P c ). This occurs at the stalling radius 

^a„=( T |^) 1/2 ^W /4 . (4) 

In a stratified medium, i? s taii is crucial in determining whether a bubble will "blow out" . In such a 
medium, the initial expansion will be spherical, as in a homogeneous atmosphere. However, if the expansion 
does not stall before proceeding beyond one scale height H of an exponential atmosphere, then it will sense 
the stratification and the top of the shell will begin to accelerate and reach an infinite height in a finite time, 
which defines blowout. Hence, the dimensionless ratio b = R stE di/H is an important diagnostic for bubbles 
in a stratified medium; we surmise that b > 1 for blowout to occur. Indeed, MacLow & McCray (1988) 
find a similar dimensionless ratio D which gauges whether or not their bubbles blow out of a finite pressure 
exponential medium. Their parameter D can be written in our terminology as 

D=(^)b* = 17.9 b*. (5) 

MacLow & McCray (1988) find that a model with D = 10 (b = 0.75) is confined but one with D = 100 
(b = 2.4) blows out of the atmosphere, consistent with our expectation that blowout occurs for b ^> 1. 

In the case of the W4 superbubble, which can be blown by stellar winds alone, the estimated wind 
luminosity of the driving cluster OC1 352 is Lq — 3 x 10 37 ergs s" 1 (NTD). Hence, for this bubble, equation 
(§) gives 

L Q \ 1/2 / n \V4 /10 4 cm- 3 K\ 3/4 /lOOpc 



b 4 \3 x 10 37 ergss-V 'l™- 3 ) \ P c /k ) \ H )' (6) 
so that any standard interstellar values for n ,P c , and H will yield b > 1. In the W4 region, probably 
b 3> 1, since the scale height H is actually much less than 100 pc (see § 2), and the ambient density no is 
actually greater than 1 cm~ 3 (see § 3). Our adopted normalizing value of P c /k also represents the high end 
of estimated pressures in the ISM (Kulkarni & Heiles 1987). 



2.2. Expansion in a Stratified Atmosphere: Application to W4 

Kompaneets (1960) developed a semianalytic solution for shock wave propagation in an exponential 
atmosphere. This solution was originally used to study blast wave propagation, but we adapt it to include 
a continuous input of energy, as also done by Schiano (1985) in his study of active galaxy winds. The 
model gives a good idea of how shock waves propagate in a stratified medium without having to solve the 
full hydrodynamic equations numerically. The shape of the bubble is determined analytically, and the time 
evolution is determined by integrating a differential equation. The spatial and temporal solutions of the 
Komapneets model for a wind-blown bubble are presented in Appendix A. In the following sections, we 
apply the solution to the W4 superbubble. 
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2.2.1. Relation of the Spatial Profiles to the W4 Superbubble 



Figure 1 shows the surface of the shock front r(z,y) (eq. |A8|) for various values of the dimensionless 
quantity y = y/H, which parametrizes the evolution of the bubble (see § A.l). The shock surface is spherical 
at early times, but becomes increasingly elongated in the vertical direction after the top surface has moved 
beyond one scale height. Several streamlines of the flow are represented by dashed lines. The elongated 
bubbles exhibit a striking resemblance to the observed W4 superbubble in the region above the star cluster. 
Figure 2 shows an H I map of NTD and Figure 3 shows the Ha map of DTS, both overlaid with the best 
fit Kompaneets model, determined as follows. The best fit is determined by finding the value of y [in the 
full range (0,2)] for which the aspect ratio matches that of the observed superbubble. The Ha map of DTS 
reveals the ionized wall of the superbubble, and shows that it reaches a maximum diameter of 3.6° at latitude 
4° and thereafter starts to close. Figure 3 shows that there is a left-right asymmetry in the structure of 
the Ha shell. The maximum radius of the shell, measured relative to the longitude of the cluster (134.7°), 
is greater towards the right (2.1°) than towards the left (1.5°); also, the shell is brighter towards the left. 
This is consistent with the interpretation that the bubble has blown into a denser medium on the left side, 
so that it has not expanded as far. When ionization bounded (see § 3.2), a shell is brighter when closer 
to the cluster. A slightly lower density toward the right-hand side may also explain why the shell appears 
to extend to somewhat greater height on that side. Since our intention is to model the mean atmosphere, 
and not such details as a left-right asymmetry, we take an average value of 1.8° for the maximum radius, 
thereby finding the Kompaneets model whose diameter matches that of the superbubble. In the same spirit, 
the top of the shell, while not clearly visible in the Ha map, is estimated to be at latitude 6.9° , based on the 
shell structure on both left and right. The star cluster is located at 0.9° latitude. Using this information, 



equations ( |A9| ) and ( A10 ), which give expressions for Z\ (22), the distance from the star cluster to the top 
(bottom) of the bubble, and r max , the maximum half-width of the bubble, are used to find the evolutionary 
stage of the bubble. Equating the theoretical and observed aspect ratios yields 

* - Mi-y/2) « 3.33. 



?"max arcsin(y/2) 1.8 

The equality is satisfied when y — 1.9837. In this calculation, we have assumed that the superbubble lies 
in a plane perpendicular to our line of sight. Our value of y is also consistent with the location of the 
bottom of the bubble, since the Kompaneets model predicts that \z 2 \/zi = 0.14 when y = 1.9837, whereas 
the observations yield |za|/^i — 0.15. The purely exponential atmosphere might be less realistic below the 
cluster, and finite external pressure might play a role there if anywhere, so it is fortunate that the model 
approximately fits the bubble here as well (see §4 for further discussion of the physics near the base of the 
bubble). 

In addition to matching the aspect ratio, the Kompaneets model also yields a value for the mean scale 
height H in the atmosphere if we know the size of the superbubble. This requires a distance, which we take 
to be 2.35 kpc, based on the distance to OC1 352 estimated by Massey et al. (1995). At this distance, the 
Ha map yields z\ ~ 246 pc and r max ~ 74 pc. Since y — 1.9837 in the Kompaneets model implies that 
z\ = 9.62 H and r max = 2.89 if, matching the Kompaneets model to the observations requires 

ff~25pc. (8) 

Even though the exact position of Z\ is is not very clear in the Ha map, the unmistakable presence of a 
maximum radius r max alone can be used to obtain a rough estimate for H, since the highly elongated nature 
of the bubble implies that r max (~ 74 pc) must be close to the limiting theoretical value ttH. Furthermore, 
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although our estimate for H is based on the Kompaneets model, it is worth pointing out that a scale height of 
this approximate magnitude is unavoidable in any model, since the current radius of the bubble (maximum 
value ~ 74 pc) must be significantly greater than H in order for the bubble to have become so elongated. 



2.2.2. The Swept-up Mass in the Shell 

The mass distribution in the shell is an important factor in determining the ionization structure of the 
atmosphere, which is calculated in § 3. The spatial profiles can be used to find the streamlines of the flow, 
since the latter are normal to the local surface at all stages of evolution (see Figure 1). Subsequently, one 
can integrate along the streamlines in the exponential atmosphere to obtain the swept-up surface (column) 
density at all points on the shell. Figure 4 shows the dimensionless surface density of the swept-up shell 
N (unit UqH) when y — 1.98 (the stage of evolution which fits the W4 superbubble) as a function of (a) 
height z/H, and (b) elevation angle 9 (in degrees) measured relative to the r-axis. The surface density of the 
Kompaneets shell is shown as a solid line, while the dashed line represents the value for expansion along a 
straight-line path to the radius s = ^/r 2 + z 2 through the same exponential atmosphere. At 9 = 0°, the radial 
(straight-line) expansion would yield a surface density N s — nor(z = 0)/3 ~ 2uqH/3, since r(z = 0) ~ 2H 
at the advanced stages of evolution of a Kompaneets bubble. The Kompaneets streamlines yield a lower N 
than radial expansion for all angles 9 > 0°, since the streamlines tend to be deflected downwards, and yield 
higher values for most angles 9 < 0° for the same reason. The crossover between the two curves occurs just 
below 9 — 0°. At 9 = 0°, the Kompaneets streamlines yield N = 0.67 UqH, versus N s = 0.69 uqH for radial 
streamlines. The Kompaneets streamlines yield values of N up to about 15% higher than the radial value 
near the bottom of the bubble, but is more than an order of magnitude below the radial value near the top. 
The continual downward deflection of streamlines means that most of the mass is not transported very far 
upwards, if at all. In the late stages of the Kompaneets model, about 55% of the matter in the shell is found 
below z = 0, about 75% is found below z = H, and nearly 90% below z = 2H. 



2.2.3. Age of the W4 Superbubble 



The time-dependent solution (§ |A.2| ) allows us to make an age estimate for any bubble if we know the 
stage of evolution (represented by y). For application to the W4 superbubble, we note that the dimensionless 
time t = 6.3 when y — 1.98. Converting to dimensional form, we obtain 



t = 1.18 



(i^)' /3 (^) 5/a (^f^)' /a ^ 



Although Lq can be estimated from direct observations of the stellar winds in the cluster stars or from the 
spectral type of the stars (NTD), and H is estimated from the size and shape of the bubble, a final estimate 
for the bubble age requires a constraint on no, which we obtain in §3 from the ionization structure of the 
atmosphere. 



3. Ionization Structure 



Along with the powerful winds which drive the W4 superbubble, the nine O stars also produce an 
extremely strong ultraviolet radiation field. The dynamics of the superbubble is dominated by the energetic 
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winds, but the optical and radio continuum appearance of the region is determined by these ionizing photons. 

The evolution of H II regions without powerful winds has been well studied (e.g., Spitzer 1978; Tenorio- 
Tagle 1982) and found to be incapable of producing extreme structures such as the W4 supcrbubble. The 
strongest constraint on dynamical H II regions is the inability to move large amounts of material at speeds 
much faster than the sound speed c s in the hot, ionized gas of temperature T <~ 10 4 K. Unlike the wind- 
blown superbubblc, where T ~ 10 6 K, an evolving H II region cannot expand more than a few tens of parsecs 
over its million year lifetime. The ionization structure of a wind-blown region is, however, produced by the 
irradiation of the swept-out cavity by the ionizing UV photons from the O stars. This allows observations 
of the ionized hydrogen to be used as a probe of the ambient atmospheric conditions. 



3.1. The Formation of H II Regions Without Powerful Winds 

Before calculating the ionization structure of the wind-blown bubble, it is relevant to review the ioniza- 
tion structure of H II regions in a variety of media. 



3.1.1. Uniform Medium 

Within a constant density environment, the ionizing photon flux from the central OB association 
produces a classical Stromgren sphere H II region with radius 

1/3 



tfst = ( T—T— 1 • (10) 



_3£* 

v 47T 7lg OtB , 

In the above, no is the number density of hydrogen nuclei inside the Stromgren sphere and olb — 2.6 x 
10~ 13 cm 3 s _1 is the recombination coefficient for hydrogen. For the nine O stars at the center of the W4 
superbubble, <!>» ~ 2.3 x 10 50 s _1 (DTS), and the Stromgren radius is 

, _ 3 n. 2/3 / , \ 1/3 

1 cm \ f $ N 



Once formed, the initial H II region expands due to the overpressure of the hot, ionized gas relative to 
the external neutral gas. 



3.1.2. Exponentially Stratified Medium 

A second useful example to consider is the instantaneous structure of the initial H II region in an 
exponentially stratified atmosphere with scale height H. Because the density drops precipitously with 
increasing z, it is possible that in the upward direction the H II region will be unbounded. If H is much 
larger than i?g t , so that the density gradient across the H II region is small, a quasi-spherical nebula is 
produced; however, if H is much smaller than R$ t , the ionizing photons are free to escape the stratified ISM 
in the upward direction. Neglecting the diffuse ionizing field, the critical ratio (Rst/ H)ait at which photons 
traveling directly upward can escape is found by equating the total number of recombinations within an 
infinitesimal solid angle d£l in the z direction with the UV flux $*: 

dtt f°° 1 

=/ nlexp(-2z/H)a B z 2 dzdn= -nla B H 3 dfl. (12) 
47T ,/n 4 
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Substituting equation ( |Io| ) into equation ( |l2] ) yields (i?st/-ff)crit = (3/4) 1 / 3 ~ 0.9. Thus, the condition for 
the initial Stromgren region to breakout is similar to that for the expansion of the wind-blown bubble, i.e., 
the characteristic radius must exceed the scale height. 

Figure 5 plots the instantaneous ionization boundary for a variety of i?st / H values showing both confined 
and breakout H II regions. In each case, the boundary of the ionized region is calculated by equating 
photoionization and recombination along each solid angle. Note the gradual elongation of the nebula in 
the upward direction with increasing Rg t /H and the shallow slope to the lower cavity when breakout is 
achieved. Like the initial Stromgren sphere in a uniform medium, these ionized regions represent the initial 
conditions for further dynamical evolution of the (overpressured) ionized gas. In extremely high-density (low 
i?St) or low scale height regions, the breakout solutions represent the initial conditions for champagne-flow 
H II regions (Tenorio-Tagle 1982), in which the ensuing dynamical evolution of the ionized gas at ~ 10 km 
s _1 is important. 

When Rst/H > 0.9, a substantial fraction of the atmosphere above the OB association may be com- 
pletely ionized and a fraction of the ionizing photons may escape to the Galactic halo. The elevation angle 
9 C above which photons escape to infinity is determined by evaluating an analogous relation to equation ( |l2] ) 
along each solid angle (see also Franco, Tenorio-Tagle, & Bodenheimer 1990). This yields 

TT 

sin(9 c = 0.9 — for R st /H > 0.9. (13) 
Rst 

When 9 C < tt/2 (Rgt/H > 0.9), the total fraction of ionizing photons that escape is 

f = \ (\-^sintf c + isin 3 # c V (14) 



Equation (14) yields limiting values / = when 9 C — n/2 (no escape of ionizing photons), and / = 1/2 when 



8 C = (the escape of all upward traveling photons) . 

In § 2, the scale height near the W4 region was estimated to be H ~ 25 pc, so that the critical value is 
R st — 22.5 pc. For the W4 superbubble <f>* = 2.3 x 10 50 s _1 , so that the critical density is n = 40cm~ 3 . 
Therefore it appears likely that during the initial illumination of the stratified ambient ISM, before any 
expansion via winds, extensive ionization would have occurred, with ionizing photons escaping to heights 
far above the cluster. In fact, we estimate no ~ 10cm~ 3 below, which gives Rst/H ~ 2.3, 8 C ~ 23°, and 
/ ~ 0.22. We emphasize that this is just the initial stage. Later, when the material is swept into a shell by 



the wind, the escape fraction / can change (see § 3.2.1) 



3.2. The Illumination of Wind-Swept Cavities 

Given the remarkable similarity in appearance between the W4 superbubble and the wind-blown bubble 
model presented in the last section, we consider the illumination of such a cavity by the ionizing photon flux 
from a central OB association. 



3.2.1. Structure of the Ionization Front 



As discussed in § 2.1, the superbubble is filled with a tenuous hot medium; this does not provide 
significant attenuation of the ultraviolet radiation field. The swept up gas from the bubble is deposited 
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along the cavity walls and is ionized by the Lyman continuum photon flux from the central OB association. 
The pressure within the superbubble is determined from the time-dependent Kompaneets solution (§ |A.2| ). 
Since the pressure within the cavity is much larger than the pressure in the ambient stratified medium (even 
when photoheatcd to 10 4 K), the swept up material, which is in approximate pressure equilibrium with the 
interior of the superbubble, is squeezed to higher density than the ambient material. The density within the 
ionized zone of the swept-up shell is calculated below. Owing to the n 2 dependence of the recombination 
rate, an important property of a wind-blown shell is that it absorbs a greater number of ionizing photons 
than the equivalent amount of matter in the undisturbed medium. The fraction of the shell which is ionized 
is easily determined by balancing ionizing photons and recombinations along a path from the central OB 
association through the shell. 

The ionization front is located where the Lyman continuum photon flux becomes negligible. We compute 
this as a function of elevation angle 8. For each 8 three possible ionization states exist: 1) The ionization 
front may occur within the high pressured swept-up shell, in which case a very dense but cold neutral 
hydrogen (possibly molecular) outer shell surrounds the ionized material. 2) The Lyman continuum flux 
may escape the swept-up shell, only to produce an ionization front further out in the undisturbed stratified 
ISM. 3) The entire column along 8 may be ionized, releasing Lyman continuum photons far beyond the 
bubble. State 2 is fairly rare in our model, because the swept-up material is very dense compared with the 
ambient exponentially stratified ISM, so that if Lyman continuum photons penetrate through the shell, they 
will likely break out of the atmosphere. 

To quantify the above discussion, consider a segment of the shell located at an angle 8 and distance 
s with respect to the central OB association. The shape of the bubble is completely specified in terms of 
the dimensionless parameter y. The surface density N(9,y) within this shell segment can be computed as 
described in § [2.2.2 and is proportional to uq H. We take the pressure within the shell to be P s (y) = 2 P(y), 



where P{y) is the average pressure in the bubble computed from the Kompaneets solution. This is a 
reasonable estimate when using a thin-shell approximation (see Bisnovatyi-Kogan & Silich 1995, p. 670). 
The number density of hydrogen nuclei in the shell is therefore 

Ps{y) h ,, 

Un = 2fc7V (15) 



Dimensional values for njj can be obtained using the dimensional pressure in the shell P s = 2P (pqLq/H ) 



4U/3 



where P has a lower bound ~ 0.02 at late times (see discussion in § A. 2 ). For example, using Tjj = 8000 K 
and no = 10 cm -3 (estimated below), we find that nu = 14.5 cm -3 at late times. The maximum thickness 
of the ionized shell perpendicular to the local tangent is 

AR = N(8,y)/n H . (16) 

For an ionizing flux $», the shell is ionized along the line-of-sight from cluster center to a depth approximated 

by 

Airnjj olb s 2 

The ionization depth perpendicular to the shell tangent, Asr, requires multiplication of equation (|l7|) by 
cos 4> where cj> is the angle between the line-of-sight direction and the surface normal, yielding 



Asr ( <fr* cos 4>kT n 



AR \2nN s P s a B s 2 



(18) 
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If Asr > AR, the entire shell is ionized and the Lyman continuum photon flux illuminates the ambient 
stratified ISM outside the bubble. If As# < AR, the ionization front occurs within the shell and a much 
denser, cold shell layer in approximate pressure equilibrium overlies the ionized layer. Note that raising the 
ionizing flux, lowering the surface density, and lowering the pressure inside the bubble enhance the ability of 
Lyman continuum photons to penetrate through the shell. Determination of the location of the ionization 
front, if it exists outside the shell, is a straightforward adaptation of the method described in § 3.1. 

The location of the ionization front determines the appearance of the superbubble in the H I observations. 
If the ionization front is always confined within the swept-up shell, then the region should appear as a closed 
bubble in the H I channel maps. However, if the Lyman continuum photons escape the shell along lines-of- 
sight above a critical 9, then the ambient medium within a cone of angle tt/2 — 9 around z will be entirely 
ionized, and the region will appear as a "chimney" in the H I channel maps. Ha observations complement 
the H I channel maps by highlighting the dense ionized gas within the shell even if the ionizing photons have 
penetrated into the low density medium beyond it. Thus, the different morphologies for the W4 superbubble 
described in NTD (H I: chimney) and DTS (Ha: closed bubble) are reconciled if the Lyman continuum flux 
breaks out toward the pole. 

The height at which breakout of the Lyman continuum photons from the dense shell can occur depends 
on the amount of material which has been swept up and on the pressure in the bubble, which sets the 
density within the ionized component of the shell. In § 2.2.% the surface density of the shell was shown 



to vary dramatically from the base of the bubble to the top, decreasing faster than the increase in surface 
area (which varies as s 2 ) because of curved streamlines. Thus, breakout is expected to occur first at the 
top. Although such quantities as the pressure and age of the bubble, and the surface density of the shell 
can be calculated in dimensionless form, where the units are the appropriate combination of no, Lq, and 
H, no dimensionless model exists when a fourth parameter is included. Therefore, specific dimensional 
models must be computed and analyzed. However, all the relevant variables except the number density are 
reasonably determined for the W4 superbubble, requiring only a one-dimensional search. The opening angle 
of the chimney observed in the H I channel maps is the basis for determination of no. Figure 6 shows the 
location of the ionization front compared with the inner edge of the cavity for various values of no- First, we 
estimate analytically the value of no at which the ionization front becomes trapped in the shell near z = 0. 
The trapping condition, Asr < AR, can be rewritten using equation ( |l8| ) as 

VKS^OLBr 

where P — 1/2 P s and cos0 ~ 1. Furthermore, on dimensional grounds one expects that 

N = Nn a H, (20) 
s = sH, (21) 
P = PipoLl/H*) 1 ' 3 . (22) 

The Kompaneets solution gives values for the dimensionless constants N, s, and P at the current epoch, 



namely N ~ 0.7 at z — 0, s ~ 2 at z — 0, and P ~ 0.02. Using these relations, equation (19) is transformed 
into the condition 



n > 
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Figure 6, which is obtained using the same normalizing values as in equation (^3j), also shows that the 
ionization front is trapped within the shell at z = when n — 5 cm -3 . As no is increased further, the 
ionization front gets trapped to progressively greater heights. For fixed no, the ionization front opens outward 
in a cone-like fashion above a critical height, as neither the shell nor the low density medium beyond it are 
able to confine the ionizing photons. Beyond a density no ~ 10 cm~ 3 , this opening point of the ionization 
front rises very slowly with increasing no, since most of the original mass is deposited within a few scale 
heights of the cluster (see the streamlines in Figure 1). We find our best fit model to the W4 superbubble 
by picking the model whose ionization front turnoff point coincides with the height of significant drop off of 
H I emission. For NTD, this occurs at about 2.9°(^ 100 pc) above the cluster (see Figure 2). This opening 
point is bracketed by the n = 10 cm -3 and n = 15 cm~ 3 curves, with no — 10 cm -3 a reasonable best fit. 
Further H I observations should be able to confirm whether the apparent drop off in H I emission continues 
at higher latitudes, which would be consistent with the idea that the ionizing photons have penetrated the 
shell at heights ;> 100 pc. It is interesting to note that the density estimate no ~ 5 cm -3 made by NTD 
by looking at H I gas adjacent to the chimney and at the latitude of the cluster is about the same as the 
minimum density allowed by our model (eq. p3|). Furthermore, it is in good agreement with our derived 
most likely value no ~ 10 cm -3 given the uncertainties in both the observations and in some of our model 
parameters. 

The fraction of Lyman continuum photons which escape the superbubble can be determined from these 
models and the analytic formulae described in § 3.1. For low density environments, little material is swept 
up by the expanding bubble, and though the bubble appears brighter than the background, the fraction 
of ultraviolet photons trapped by the shell is insignificant; the escaping fraction may be deduced from the 
arguments for a static exponential atmosphere, i.e., equations (|l3| ) and (fl4|). However, when the density of 
the medium is large enough, the exponential atmosphere becomes insignificant for this calculation and the 
fraction of photons that escape is determined by the absorption within the illuminated shell. For our best 
fit model (no ~ 10 cm -3 ), the fraction of solid angle through which ionizing photons escape (when y — 1.98) 
is approximately 15%. At these opening angles, the shell is at a relatively high latitude and contains little 
mass, so that nearly all the ionizing flux is transmitted. Therefore, we calculate that a flux of ionizing 
photons ~ 0.15 $* ~ 3 x 10 49 s" 1 from the cluster OC1 352 escapes to higher latitudes. 

Below z ~ 100 pc, and particularly in the downward direction [z < 0), this model predicts that the 
stellar cluster does not completely ionize the dense wind-swept shell. There should then be a cooler, even 
more dense shell outside the ionization front, again in approximate pressure equilibrium. It is therefore 
interesting that the IRAS observations show peak emission just outside the ionized shell along the lower 
part of the cavity (Heyer & Terebey 1998). For maximum conversion of stellar into reradiated infrared 
luminosity by dust, the shell should have substantial optical depth. This can be checked by noting that the 
column density at z ~ Z2 is ~ no H ~ 10 21 cm -2 (see Figure 4) for our standard values of no (— 10 cm -3 ) 
and H (~ 25 pc). For a normal dust to gas ratio, this implies that Ay ~ 0.5 mag, yielding significant 
self-obscuration of our view into the lower bubble. Accompanying this optical depth, we have the possibility 
of molecular gas (which might be swept up, having survived the earlier evolution, or have reformed). Heyer 
& Terebey (1998) also show CO emission confined along a thin layer adjacent to the bottom edge of the 
superbubble, indicating the presence of dense swept up molecular gas. 
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3.2.2. Emission Measure 



Considering the evolutionary scenario for a wind-blown superbubble outlined in § 3.2.1, it is straight- 
forward to produce models of the emission measure profile, J n 2 dl, where the integral is taken through the 
bubble along our line-of-sight, for a variety of no values (Figure 7a-c). In each panel six epochs are shown 
corresponding to y = 0.5 — 1.98 and the emission measure is presented as a grey scale, with darker regions 
corresponding to higher emission measure. A cavity in H I channel maps should occur wherever emission 
is observed; in these regions the gas is ionized. Each figure also shows the instantaneous location of the 
ionization front separating the ionized from the neutral gas. Thus, when no in the ambient ISM is low, as 
in Figure 7a, the region always appears as a giant cavity in H I . The ionization front occurs in the ambient 
medium outside the shell, so that in Ha or radio continuum the region shows a thin-shelled bubble with a 
very bright halo beneath. As no is increased, the H I chimney narrows and the bright Ha halo disappears 
from beneath the bubble. In the highest density case presented, (no — 10 cm -3 ), the upper shell is opaque to 
the ionizing photons for a time, once sufficient mass is swept up into it, but eventually becomes transparent 
when the internal pressure (and the density of the shell) has dropped sufficiently. 

The Ha shell thickness observed by DTS falls in the range 0.2° to 0.5° at all heights. Our model predicts 
a thick shell near the bottom of the bubble, and an extremely thin shell (much thinner than observed by 
DTS) in the upper portion of the bubble (Figures 7a-c). The former is due to the large amount of mass 
swept up in the increasing exponential atmosphere below b = 0.9° (z = 0), which is probably not a good 
representation of the atmosphere there, and the latter is due to the uniform pressure assumed in the bubble 



(see § A. 5 for further discussion). A much lower pressure should exist near the top of the bubble at this 
stage, allowing for a thicker shell there. However, our model should give a reasonable estimate for the shell 
thickness near z = 0. The maximum shell thickness there can be found by combining equations (|l^), jl^), 
©, and (H), yielding 

AR = N n 2 n /3 H 7 / s kT 



2/3 



A>\/0.02\/ no \ 2 /3/ H x 7/3 / 9 - 1 n37 — N 2/3 



= 1L1 ^J^Jlio^J fej { Z-T^J pc - ™ 

At a distance d ~ 2.35 kpc, the standard value obtained above yields 

A R 

A9 = =f ~ 0.27°, (25) 
a 

which falls within the range observed by DTS. 

A direct comparison between the Ha emission produced by the best-fit model of the ionized wind- 
swept bubble (no = 10 cm -3 ; last frame of Figure 7c) and the Ha observation of DTS reveals the merits 
of the simulation. Conversion between emission measure and rayleigh flux units (R) is straightforward: at 
Tn = 8000 K, 1 R = 2 pc cm -6 . DTS compute an extinction at Ha of 2.1 mag toward the star cluster 
which we apply as a uniform obscuring medium to the model bubble. Also, as shown in Figure 8, the 
intensity profile of the DTS observation (averaged over 9 pixels or 9 pc) along the symmetry axis r — 
reveals a high background flux below the bubble and a low background flux above the bubble (solid line). 
We empirically model this background as an exponential distribution with a characteristic scale of 3.3° and 
a 140 R normalization at b = 0.9° (the latitude of the stellar cluster) before attenuation (dashed line in 
Figure 8). The best-fit model including the empirical background and attenuation to the bubble is also 
shown in Figure 8 (thick line). The overall profiles match remarkably well despite the many internal knots 
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of emission appearing in the observations which have not been modeled. The brightness of the model at 
the bottom of the shell (z <J —20 pc) is significantly greater than the observed Ha flux. This deficit of 
observed Ha flux was also noted by DTS, and arises in comparison with any model which assumes that 
the bottom of the shell is ionization bounded, since the Ha brightness is then determined by the ionizing 
flux. A possible explanation for the discrepancy is extinction due to dense swept-up neutral gas outside the 
ionization front (Ay ~ 0.5 mag; see earlier discussion). However, the edge-brightening enhancement at the 
bottom shell appears reasonable as does the slow decline in brightness with height. The intensity declines in 
two stages, separated at z ~ 100 pc, which marks the position within the model where the ionization front 
breaks through the wind-swept shell. Both components appear visible in the Ha observations. 

In Figure 9, an intensity profile is produced across the bubble at a height z = 50 pc. The model has 
been attenuated in the same manner as discussed for Figure 8 (at this height the background flux is 100 
R in the empirical model). The Ha data reveals enhanced emission within the bubble comparable to the 
level expected from the model, but the two wall features are not observed. An intensity peak of the correct 
magnitude is observed at r — —40 pc, and is consistent with the fact that the Ha shell is pinched at low 
latitudes (Figure 3). As noted above, the thickness of the observed peak is greater than derived in the 
model, indicating that at this latitude the pressure in the bubble may not be as high as expected from a 
simple Kompaneets solution. The bubble edge is not apparent on the right hand side; however, the intensity 
does decrease significantly where the edge should be. Part of the difficulty in matching the edge-brightened 
intensities may be the requirement, satisfied in the model but perhaps not in reality, of significant coherence 
along our linc-of-sight. If the thin bubble wall has non-axisymmetric structure, then the intense edge emission 
will be spread over a larger area and the bright cusp will not be observed. 



4. Discussion 

Given the best fit parameters no ~ 10 cm" 3 , Lq ~ 3 x 10 37 ergs s _1 , and H ~ 25 pc for the W4 
superbubble, equation (||) yields the best-fit age estimate 

t ~ 2.5 Myr. (26) 

This is less than the estimated main-sequence lifetime (3.7 - 4.3 Myr) of the O stars in the cluster (NTD), 
and is therefore consistent with NTD's interpretation that the superbubble is blown by the combined winds 
of the young cluster in which no supcrnovac have yet occurred. The cluster age has also been estimated by 
various authors to be < 2.5 Myr (see DTS). The superbubble is on its way to blowing out of the Galactic disk 
unless the atmospheric structure changes significantly at greater heights. Given our best fit parameters, the 
blowout parameter (defined in § |]l]) b ~ 25 (see eq. 0j), so that use of the Kompaneets model is justified. 
A quick estimate of the superbubble age can also be made using an important feature of expansion in a 
stratified atmosphere pointed out in § A. 2; the base radius at z = 0, about 50 pc in this case, can be equated 
to R s in the spherical solution (eq. [jl|), and our best estimates for no and Lq used to obtain t ~ 2.1 Myr. 
This is in reasonable agreement with the value (eq. J2(|) determined directly from the Kompaneets model. 

The shape of the W4 superbubble provides an important surprise: the ambient interstellar gas around 
the cluster OC1 352 was very strongly stratified in the vertical direction. The estimated scale height H ~ 25 
pc is based on the best fit Kompaneets profile and the estimated distance 2.35 kpc to the cluster. This is to be 
compared with typical estimates H ;> 100 pc for interstellar gas near the Galactic plane based on emission 
from H I or CO. For example, the FCRAO CO survey of the outer Galaxy (Heyer et al. 1998), where 
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the Perseus arm lies, finds an overall mean CO scale height of 113 pc[J. Our scale height determination 
is, in contrast, based on a dynamical model of an observed superbubble, and is applicable only to the 
immediate vicinity of the stellar cluster OC1 352 and the W4 H II region. Our scale height estimate would 
increase linearly in proportion to the cluster distance if it were actually farther than the 2.35 kpc recently 
determined by the spectroscopic parallax method (Massey et al. 1995). However, we do not expect the 
actual distance to be much different, since many previous estimates for the distance to OC1 352 or the 
W3/W4/W5 H II regions have, using a variety of techniques, yielded values in the range 1.6 to 2.7 kpc (see 
NTD). Similarly, the estimate H ~ 25 pc cannot be explained away due to any of the simplifications of the 
Kompaneets model. All models for bubble expansion, whether analytic, semianalytic, or numerical (e.g., 
Kompaneets 1960; MacLow & McCray 1988; Tomisaka & Ikeuchi 1986; MacLow et al. 1989), show that 
the bubble remains nearly spherical as long as the ratio of the radius in the plane z = to the scale height 
R/H < 1. Significant elongation takes place only when R/H ;> 2, i.e., the bubble has had a chance to fully 
sample the stratified atmosphere. The W4 superbubble is currently highly elongated and has a radius in the 
plane R ~ 50 pc, hence one can conclude that H < 50 pc, and most likely H ss 25 pc, using very simple 
physical arguments. 

The relatively small value of H in comparison to the mean ISM value (> 100 pc) implies that significant 
vertical compression of the interstellar gas has taken place here, perhaps not such a surprising finding for 
a star-forming region. Consistent with this, the most likely value for no (~ 10 cm -3 ) is also significantly 
greater than the mean ISM value (~ 1 cm -3 ). It is interesting to note that the total surface density of gas, 
which is proportional to uqH, is about a factor of two greater than implied by the mean ISM values. This 
is consistent with the idea that some lateral inflow of gas has also taken place during the assembling of the 
molecular cloud. 

The Kompaneets model idealizes the atmosphere as being purely exponential, so it is all the more striking 
that the y — 1.98 profile provides such a good approximation to the overall shape of the superbubble (Figure 
3). Thus, the exponential atmosphere and scale height H ~ 25 pc should be reasonable approximations 
for the mean atmosphere into which the superbubble has expanded. At large heights, towards the top of 
the superbubble, it is possible that the stratification is not so great. Although there is relatively little mass 
there, it may be more than implied by an exponential atmosphere with H ~ 25 pc. Possible deviations from 
the assumed atmosphere, e.g., the extent to which the medium with H ~ 25 pc can be embedded in one with 
larger H at high latitudes and/or adjacent longitudes, can be best addressed with future modeling using the 
thin shell approximation or numerical models. However, we note that various models of bubble expansion 
(e.g., MacLow & McCray 1988) have shown that a bubble changes shape dramatically when it travels from 
a relatively low scale height atmosphere (e.g., the H I cloud layer) to one with greater scale height (e.g., the 
H I intercloud layer); that is, the bubble radius expands to become comparable to the local scale height. 
This sort of expansion has not yet happened to the W4 superbubble. Hence, we conclude that the effect of 
extra mass near the top of the bubble cannot be very great. 

Below the star cluster, the density probably decreases eventually, certainly on passage past the midplane 
of the Galaxy, unlike the increasing exponential function in the Kompaneets model. However, the presence 
of a decreasing density gradient in the downward direction must be much less pronounced than that above 
the cluster, or else the superbubble would blow out in both directions. The finite pressure of the external 



1 One should note that this is based on kinematic distances determined by the standard circular rotation curve of the Galaxy, 
which can overestimate the distance for gas in a spiral arm, thereby also overestimating the scale height. Heyer & Terebey 
(1998) point out this effect; the kinematic distance to OC1 352, based on associated molecular gas, ranges from 3.5 to 4.0 kpc, 
while the spectroscopic distance to OC1 352 is only 2.35 kpc (Massey et al. 1995). 
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medium may also play a role in confining the downward expansion. 

Deviations from a pure exponential atmosphere are already evident in Figures 2 and 3. Figure 3 reveals 
that the superbubble is narrowest at a latitude almost equal to that of the star cluster. Such a pinch is 
usually caused by a local maximum in the density (e.g., MacLow & McCray 1988). There is clear evidence 
for high density regions which inhibit the expansion immediately adjacent to the cluster; CO and IRAS 
images reveal a dense cloud near the left (east) pinch and the right hand (west) side is certainly affected by 
the presence of the W3 molecular cloud. 

Further evidence for a more complicated atmosphere is provided by the H I and CO maps (NTD; Heyer 
et al. 1996), which reveal embedded clouds in the cavity. These clouds are most likely remnants of an 
initially clumpy medium; a superbubble can move around dense clumps and continue to expand (MacLow 
& McCray 1988). The CO maps reveal two cometary shaped clouds (Heyer et al. 1996), and NTD find a 
V-shaped filament in H I (visible in Figure 2), with the higher latitude CO cloud at its base. The V-shaped 
"streamers" are suggestive of material being stripped from the molecular cloud by a free-streaming wind, 
but Heyer et al. (1996) argue that the cometary shape of the CO clouds is likely due to the sculpting 
effects of the ionizing UV radiation field of the O stars. Our model supports the view that the CO clouds 
are embedded in a hot rarefied gas in which the ionizing UV radiation from the central source propagates 



freely, i.e., region (b) in the description of a wind-blown bubble given in § 2A. The region (a) consisting 
of free-streaming wind remains relatively small during the evolution due to the high pressure in the hot 
shocked region (b). We estimate the current radius of this region to be Ri ~ 10 pc (see Appendix B). For 
comparison, the bottom of the lower latitude embedded CO cloud observed by Heyer et al. (1996) is about 
16.5 pc above the cluster. We note that photoevaporation of these clouds will act to cool the shocked gas in 
the bubble interior (McKee, Van Buren, & Lazareff 1984) and slow down the expansion of the bubble. This 
will tend to increase the inferred age of the bubble given by equations (m and 



The thickness of the ionized shell in our model agrees well with observations of the superbubble near 



0, as shown in § 3.2.2. However, the shell thickness in our model becomes progressively thinner at 



high latitudes than revealed by the Ha map of DTS. This can be due to several reasons. As explained in 



§ A. 2, our estimate of the pressure applies only near the base of the bubble. In reality, a pressure gradient 
must be present in the bubble at this stage, so that the pressure is much lower near the top, allowing for 
a thicker (but less bright) shell there. Additionally, a swept up magnetic field may contribute significant 
pressure within the shell, and prevent it from being compressed to the level in our model. A tangential 
magnetic field within the shell might also explain another puzzling feature of the observed superbubble. The 
Kompaneets model, and indeed all expansion models, predict that a bubble is in the rapidly accelerating 
blowout phase by the time it is so elongated (recall the decelerating expansion only lasts while the bubble 
is within one or two scale heights above the plane, and still has a nearly spherical shape). Hydrodynamic 
simulations of expanding superbubbles show that the shell fragments during this stage due to a Raylcigh- 
Taylor instability (MacLow et al. 1989). Yet, the Ha map of DTS shows that the W4 supershell is relatively 
smooth. A plausible explanation for the suppression of the Rayleigh- Taylor instability is the stabilizing effect 
of a swept-up tangential magnetic field in the shell (e.g., the simulations of Stone & Norman 1992). We 
investigate this issue further in a future paper; preliminary results can be found in Komljenovic, Basu, & 
Johnstone (1999). 

While magnetic fields might explain some features of the W4 superbubble, we caution that they can also 
raise complications. For example, a dynamically significant magnetic field aligned predominantly parallel to 
the Galactic plane would prevent the bubble from attaining its present elongated shape, both by potentially 
increasing the scale height H and also by exerting a force that inhibits expansion perpendicular to the plane. 
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Indeed, simulations of bubble expansion in a stratified medium with a dynamically significant field parallel 
to the Galactic plane show that the bubble expands preferentially along the Galactic plane (Tomisaka 1992), 
unlike the W4 superbubble. 

Our model for the ionization structure of the bubble also predicts that a significant fraction (~ 15%) of 
the ionizing UV photons from the cluster escape above the bubble. Although the top of the W4 superbubble 
is at a height z\ ~ 246 pc, and we have no information about the density structure above the superbubble, 
it is quite likely that if the ionizing photons can penetrate the dense upper shell of the bubble, they will 
then continue to ionize the higher regions of lower density. The escape of ionizing photons is ultimately 
traced to the relatively low scale height H of the local medium in comparison to the initial Stromgren sphere 
radius of the combined O stars. If relatively low values of H, on the order of the value found in this paper, 
are common around major star-forming regions and young star clusters, then the consequent penetration of 
ionizing UV photons far above the Galactic plane may be a natural explanation for the extensive presence 
and large scale height (> 1 kpc) of free electrons in our Galaxy (e.g., Reynolds 1989, 1991). 

5. Summary 

We have used a model for the dynamical expansion and ionization structure of a superbubble to gain 
insight into ISM conditions near the young cluster OC1 352. The Kompaneets solution provides a reasonable 
fit to the observed shape of the H I cavity and Ha shell of the W4 superbubble, and when combined with a 
distance to the cluster, yields an estimate for the mean scale height H of the local atmosphere. Remarkably, 
we find that H ~ 25 pc in this star-forming region, significantly less than the mean ISM value near the 
Galactic plane. This relatively low value is an unavoidable conclusion given the significant elongation of 
the superbubble. The Kompaneets solution also allows us to estimate the current pressure near the base of 
the superbubble. This estimate, along with properties of the ionized shell of the superbubble, enable us to 
constrain the value of the original ambient density no near the cluster. The total wind luminosity of the 
massive stars L can be estimated from observations and/or the spectral type of the stars (Normandeau et 
al. 1996). Therefore, we can use the best values for no (~ 10 cm~ 3 ), Lq (~ 3 x 10 37 ergs s _1 ), and H (~ 25 
pc) to find a dimensional value for the age of the superbubble. Our estimated age of 2.5 Myr is consistent 
with the age of the cluster and the idea that the superbubble is blown by the combined winds of the cluster 
stars, among which no supernovae have yet occurred. 

We have calculated the ionization structure of a stratified medium, both before and during the expansion 
of a superbubble within it. The combined photoionizing ability of the nine O-type stars is sufficiently high 
and the ambient scale height H of the medium sufficiently low that initially the ionizing photons escape to 
heights far above the cluster. The expansion of the superbubble creates a dense shell, which could even trap 
all the ionizing photons for a time. However, during the late stages, when the bubble pressure (and therefore 
the shell density) has dropped sufficiently, the ionizing photons penetrate beyond the upper portions of the 
shell (of lower mass surface density) and escape to greater heights; the shell still traps the ionizing photons 
near the base (z <^ 0). This means that while the upper shell will be an enhanced region of Ha emission, 
the superbubble will have no visible lid in H I . Thus, while the apparent closure of the Ha shell (Dennison 
et al. 1997) places a damper on the Galactic chimney hypothesis (at least at the present epoch), the W4 
superbubble may continue to appear as a chimney when observed in H I at higher latitudes than in the 
CGPS Pilot project (Normandeau et al. 1996). Finally, if the presence of a relatively low scale height of gas 
is common near an OB cluster, the consequent escape of ionizing photons may provide a natural explanation 
for the extensive presence of ionized gas in our Galaxy. 
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APPENDICES 



A. KOMPANEETS MODEL 
A.l. Spatial Solution 

The Kompaneets model assumes that the energy source is located at z — in an exponential atmosphere 
(here the ISM) with density 

p{z) = po exp(-z/H), (Al) 

where H is the scale height. Note that z = is not necessarily the Galactic midplane. The model assumes 
that the internal pressure of the bubble dominates any external pressure: the atmosphere is pressure-free. 
This means that all bubbles will continue to expand in the z direction; furthermore, the top of the shell will 
reach infinite z in a finite time, i.e., undergo blowout. Physically, we can use the Kompaneets model when 
we are confident that the external pressure P e will not confine the bubble before the upper shell begins to 
accelerate, i.e., the parameter b (eq. ||) must be significantly greater than one. This condition is certainly 



satisfied by the W4 superbubble, as shown in § 2.1 



The internal pressure of the bubble is taken to be constant, and is 

P=(7-l)^p (A2) 

where 7 is the ratio of specific heats, E t h is the thermal energy of the bubble, and the volume f2 of the 
remnant is defined in cylindrical symmetric coordinates (r, z) by the integral 

rzi 

n = ir r 2 ( Zl t)dz, (A3) 

J Z 2 

in which z\ and z^ (< 0) are the top and bottom of the remnant, respectively. The extremely high temperature 
of the shocked gas, and hence short sound crossing time, is the reason that the internal pressure is assumed to 
have readjusted to an isobaric state at all times. Since the internal pressure dominates the external pressure, 
the expansion speed is that given by the Hugoniot conditions for a strong shock, 



The expansion is assumed to be directed normal to the local surface of the remnant. In the Kompaneets 
approximation, the mass swept up by the outer shock is assumed to reside in a thin shell behind the front, 
although the effect of its inertia is not included in the model. 



Using equation (A4) and the assumption that the shock front moves normal to itself everywhere, Kom- 



paneets (1960) derived an equation for the evolution of the shock front (see also Bisnovatyi-Kogan & Silich 
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1995): 



Or 
dy 



Po 



0. 



In equation (|Aq), y is a transformed variable (with units of length) defined by 




The thermal energy Eth of the bubble is calculated from the differential equation 



dE th 
dt 



dfl 



(A5) 



(A6) 



(A7) 



This equation differs from that used in the blast wave formulation of the problem, in which E'tii is a constant 



fraction of the energy deposited in an initial blast. Equation (A7) assumes that the wind luminosity is 
thcrmalized at the inner shock front R\ , and that the only energy loss is due to work done against the thin 
shell. The effect of radiative and evaporative cooling of the bubble interior can also be added here if desired. 



Kompaneets (1960) showed that equation (A5) could be solved analytically by separation of variables, 
yielding the relation 

- l / .2 

,-z/H 



r(z,y) = 2H arccos 



2 \ 4ff 2 



(A8) 



Equation (A8) describes a sequence of shapes for the shock front, which change as the parameter y varies 
from to 2H , at which time the top of the remnant formally reaches infinity. Note that one can discuss the 
structure of the shock front without actually solving for the explicit time-dependence y versus t (§ |A.2|) from 



equations ( A6) and ( A7). The top and bottom of the remnant, where r = 0, are located at 



z 1 , 2 = -2H MlT^:). 



The maximum radius of the bubble, where dr/dz = 0, is 



2H arcsin ( ) 

\2HJ 



(A9) 



(A10) 



The above equations show that by the time of blowout (y — 2H), the bubble has radius in the plane 
r(z = 0) = 2.09-ff, maximum radius r max = nH, and has penetrated downward into the atmosphere of 
exponentially increasing density to a location zi = — 1.39-ff. 



Equation (A8) illustrates an important property of the Kompaneets model; the spatial solution is 
independent of the time evolution. The shock front evolves through the sequence of shapes given by equation 



(AS) which is a consequence only of the atmospheric structure. However, the rate at which it evolves does 
depend on the details of the energy input into the bubble, i.e., whether the energy is input at one instant or 
in a continuous fashion, and at what rate. Hence, the shape of an observed bubble can supply information 
about the ambient atmosphere independent of any knowledge of the energetics of the driving source. 



A.2. Time Evolution 



A numerical integration of equations ( [A6] ) and (A7) yield y(t) and Eth(t), which implicitly give the 
time evolution of the bubble. The solution for a wind-blown bubble in an exponential atmosphere depends 
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on only three parameters: the scale height H, the density near the source pq, and the wind luminosity 
Lq. Hence, the problem is solved most naturally in a system of units which are formed from these three 
parameters. The unit of length is then H, the unit of mass is PqH 3 , and the unit of time is (pqH 5 / Lq) 1 / 3 . 
The units for various physical quantities under this system are listed in Table 1 . We obtain the dimcnsionlcss 



solution by integrating the dimensionless form of equations flAq ) and (A7). The dimensional solution for 
any combination of pq, Lq, and H can be obtained from this single solution using the units shown in Table 
1. In the following discussion, dimensionless variables are overlaid with a tilde sign. Figure 10a shows the 
parameter y versus time t in a solid line. The integration is carried up to y = 1.98, the stage at which the 
spatial profile matches the W4 superbubble. We do not integrate all the way to the limit y — 2 due to the 
artificial nature of the solution at that time; the top of the remnant reaches infinity at infinite speed and 
the pressure drops to zero. For comparison, the dimensionless radius of the spherical expansion solution 
R s = 0.76 1 3 / 5 (from eq. Q) is plotted alongside in a dashed line. The relationship between y and R s helps 
to assign some physical meaning to y. For early times, when the expansion is nearly spherical, inserting the 



spherical homogeneous solutions for R s (eq. 0]) and E t h (eq. into equation (A6) yields 

y = 0.77P /5 (All) 

if we take 7 = 5/3. Therefore, y is essentially equal to R s at early times. However, Figure 10a also shows 
that y ~ R s throughout most of the evolution. Although y begins to fall below R s when the bubble has 
expanded beyond one scale height, it is only 14% below i? s when y — 1.98. Figure 10b shows the distance 
traveled by the bubble shell in three directions: 1) upward along the z-axis, i.e., z\, 2) downward along the 
z-axis, i.e., — 22, and along the r-axis, i.e., r(z — 0). The dimensionless spherical solution R s is again given 
for comparison in a dashed line. The top of the bubble z\ eventually moves much farther than R s since it 
moves into a progressively lower density environment, while the bottom moves less far than R s as it moves 
into a progressively higher density environment. With constant density along z = 0, the time evolution of 
r(z — 0) closely follows that of R s , even though the curved streamlines mean that r(z = 0) does not follow a 
mass element in a Lagrangian fashion. There is an even closer correspondence between y and r[z = 0); even 
at the end of the integration, the two values differ by only 4%. Altogether, we find the important relation 

y~r(z = 0,y)~R s . (A12) 

Therefore, the measurement of the radius r(z = 0) of any observed bubble can yield quick estimates for y 
and even the age t (if one has estimates for Lq and po) since R s can be converted to a time through equation 
(i- 

The evolution of the thermal energy E t h(t) is shown in Figure 10c. For comparison, the dimensionless 
spherical solution (from eq. [0) is shown in dashed line. Interestingly, £th(i) remains very close to the linear 
time-dependence of the spherical solution even after there is significant elongation of the bubble. This means 



that the term Pdfl/dt in equation ( A7) remains nearly constant even when the bubble is becoming elongated. 
At late times, E t h(t) does fall significantly below the spherical expansion value and even begins to decrease. 
This decrease is due to the rapid expansion of the bubble in the late phases; dVL/dt increases without bound, 



approaching an infinite value at y — 2, causing the second term in equation ( A7) to dominate the first, making 
dEth/dt negative. However, the Kompaneets approximation breaks down when the expansion speed exceeds 
the internal sound speed, so we cannot trust that the decrease of -Eth(^) will occur in a more realistic model 
for this stage. Figure lOd shows the pressure P versus t. Again, the dashed line shows the dimensionless 
spherical solution from equation (|^). As with other variables, the pressure follows the spherical solution for 
much of the evolution, but it does drop significantly below it during the late stages. This stage coincides 
with the rapid expansion of the bubble, when one of the key assumptions of the Kompaneets approximation, 
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uniform pressure within the bubble, begins to break down. Uniform pressure certainly cannot be maintained 
when the expansion speed exceeds the internal sound speed. A pressure gradient is then set up within the 
bubble. Hence, the pressure at the base of the bubble, near the cluster, will not decrease to the values shown 
in Figure lOd. 

We estimate the current pressure at the base by equating it to the mean pressure P at a time when 
rapid acceleration begins and the expansion speed approaches the internal sound speed. Beyond this critical 
time, we assume that the pressure at the base will not be able to readjust to a new value in the short time 
available. This gives us an estimate of the pressure at the base without going through the complexities of a 
full numerical solution of the hydrodynamic equations. We estimate the time at which the base pressure is 
frozen at the mean value by comparing dz\ /dt with the internal sound speed c s ,i = \J 'kT^ / /itoh- The bubble 
temperature Tb is estimated in Appendix B. Figures lOe and lOf show the speed of the top boundary dz\jdt 
versus t and z% respectively; the latter gives an idea of where, rather than when, the acceleration phase 
occurs. Overplotted on both curves, in dashed lines, is the internal sound speed c Sj i = c Sj i/[w], where [v] is 
the unit of velocity given in Table 1. The standard values Lq ~ 3 x 10 37 ergs s _1 and H ~ 25 pc are used to 
evaluate [v], and no is taken to be 1 cm" 3 . We use dimensional values for Lq, H, and no in this case since 
the temperature depends implicitly on the thermal conductivity C within the bubble. Hence, the sound 
speed depends on a fourth independent parameter and does not have a unique dimensionlcss value. Figures 
lOe and lOf show that the Kompaneets approximation breaks down after the top of the bubble departs from 
the i -2 / 5 deceleration of the spherical solution and goes into the acceleration phase. The crossing of the 
two curves takes place at dimensionless time t — 4.9 for Uq = 1 cm -3 , when the dimensionless pressure is 
P ~ 0.02. Although this crossing would occur at a slightly later time, and lower pressure, if we use a higher 
value for no, we use this dimensionless critical pressure in our analysis, since the combination of near-sonic 
velocity and very rapid acceleration make pressure readjustment unlikely after this time even for somewhat 
higher no in the range 1 to 10 cm -3 . In § 3, we use the lower limit on the base pressure, P ~ 0.02, as a 
parameter in our model for the ionization structure of the atmosphere. 



B. SOME INTERIOR PROPERTIES OF THE W4 SUPERBUBBLE 

B.l. The Interior Temperature 

The bubble temperature Tb is obtained from the model of Weaver et al. (1977). Their equation (37) 
for the interior temperature T in the bubble can be rewritten as 



T = 63.6 {'-f) ' (1 - r/Rs) 1 ^ K. (Bl) 



It is valid for r < R s , where R s is given by equation (|lj). Equation (Bl) is obtained under the assumption 
of balance between the outward thermal conduction flux and the inward mechanical energy flux due to 
mass evaporation from the cool dense shell (see Weaver et al. 1977). Though our bubble is highly elongated, 



equation (Bl) should give a reasonable estimate for the temperature near z — since we find that the solution 
there is in many respects similar to the spherical solution. Normalizing to the standard wind luminosity and 
current radius R s ~ r(z — 0) ~ 50 pc of the W4 superbubble, we find a central temperature (which is a 



representative value through most of a spherical bubble due to the weak radial dependence of eq. |B1|) 

2/ 7 /Kn™\ 2 /7 



Tl - 5xl0 6 ( Lo \ / 50 pc 

Jb ~ 5Xl ° UxH^ergss-iJ ^ R, 



K. (B2) 



-21 - 



B.2. The Inner Shock Radius 

We estimate the inner shock front radius Pi from the ram pressure of the free wind 

p ~=**=^*--U-> (B3) 

where p w is the wind density, v w is the terminal velocity of the wind, M is the mass outflow rate, and the 
luminosity L — l/2Mv^. Equating this to the pressure Pb inside the region of hot shocked gas yields a 
radius 



{2L Mf/ i 
(47rPb) 

Now, the pressure in the bubble is written in dimensional form as 



* = (B4) 



P h = P{p Ll/H i ) 1 '\ (B5) 



where P is determined from the time evolution of the Kompanccts solution. In § A. 2 we estimated that 
P ~ 0.02 at the current epoch, hence using this as well as the most likely values uq ~ 10 cm~ 3 , Lq ~ 3 x 10 37 
ergs s _1 , and PT ~ 25 pc for the W4 superbubble, we find that 

P h /k ~ 1.1 x 10 5 cm" 3 K, (B6) 



where A; is the Boltzmann constant. Using this value in equation (BJ) with a total mass loss rate for the 
cluster stars M ~ 10 -5 Mq yr" 1 (see NTD), we find 

Pi ~ 10 pc. (B7) 



Hence, the hot shocked gas occupies nearly the entire volume of the superbubble, and since Pi < H, we are 
in retrospect justified in using spherical geometry to calculate Pi. 
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Fig. 1. — Evolution of the shock front in the Kompanects solution, as the bubble expands in an exponential 
atmosphere p = poexp(— z/H) and eventually blows out. Solid lines show the position of the shock front 
at seven successive times characterized by increasing values of the dimensionless parameter y ( = 0.5, 1.0, 
1.4, 1.7, 1.9, 1.98, and 2.0). Dashed lines show the streamlines of the flow for various points along the shock 
front. 
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Fig. 2.— H I channel map showing the location of neutral hydrogen in the vicinity of W4 and revealing the 
empty chimney discovered by NTD. Lighter shadings correspond to brighter H I emission. Overlaid on the 
figure is a Kompanects solution with y = 1.98 and H = 25 pc, centered on the star cluster OC1 352. 
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Fig. 3. — Ha map showing the location of ionized hydrogen in the vicinity of W4 and revealing the ex- 
tended superbubble discovered by DTS. Lighter shadings correspond to brighter Ha emission. Overlaid is a 
Kompaneets solution with y = 1.98 and H — 25 pc, centered on the star cluster OC1 352, as in Figure 2. 
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Fig. 4. — Surface density N (unit uqH) in the swept up shell, when y = 1.98, versus (a) height z/H, and 
(b) angle 0, in degrees, measured relative to the r-axis (z = 0). The dashed line shows the surface density if 
the expansion had occurred along straight lines. 




Fig. 5. — Location of the initial ionization front around an H II region formed within an exponential 
atmosphere. The computed models, from the innermost curve, are for Rst/H = 0.1, 0.3, 0.5, 0.7, 0.9, 1.0, 
2, 3, 4. 
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Fig. 6. — Location of the ionization front around a wind-swept bubble assuming the Kompanects solution 
(y = 1.98), the stellar properties outlined in the text, and varying only no, the number density in the ionized 
exponential atmosphere (H — 25 pc) at z = 0. From bottom to top, no = 1, 5, 10, 15, 20 cm~ 3 . The dashed 
line plots the inner boundary of the ionized wind-swept shell. 




Fig. 7. — (a). Caption on p. 33. 




Fig. 7.— (b). Caption on p.33. 




Fig. 7.— (c). Caption on p. 33. 
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Fig. 7. — Relative emission measure across an expanding bubble ionized by an internal cluster of stars with 
properties outlined in the text. Within an exponential atmosphere (H = 25 pc) six stages in evolution are 
represented (y = 0.5, 1.0, 1.4, 1.7, 1.9, 1.98). Distances are labeled in pc on both axes, with the star cluster 
located at (0,0). (a) n = 1cm -3 . Note that the ionization front is always well separated from the bubble, 
(b) no = 5 cm -3 . In this case the increasing column density at the base of the bubble confines the ionization 
front at late stages, (c) no = 10cm~ 3 , representing the best case scenario for the W4 superbubble. Note 
that the ionization front is initially open towards the top, is confined by the high density swept-up shell at 
intermediate epochs, and then breaks out when the bubble gets elongated. 




log I(R) 

Fig. 8. — Ha intensity profile along the symmetry axis (r = 0) of the superbubble. Thick line: model 
described in Figure 7(c) after conversion to raylcighs, addition of background flux (dashed line) as described 
in text, and extinction by 2.1 mag. Thin line: as observed by DTS along vertical cut centered on the cluster 
OC1 352. 
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Fig. 9. — Ha intensity profile across the model supcrbubblc at a height z = 50 pc above the cluster. Thick 
line: model described in Figure 7(c) after conversion to rayleighs, addition of background flux (see text), 
and extinction by 2.1 mag. Thin line: as observed by DTS. 
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Fig. 10. — Caption on next page. 
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Fig. 10. — Time evolution of various dimcnsionless variables in the wind-blown Kompaneets solution (vari- 
ables in (a) - (e) plotted versus the dimensionlcss time t). (a) y. For comparison, the dimensionless solution 
for the radius of a spherical bubble in a uniform atmosphere, R s , is shown with a dashed line, (b) Distance 
from the wind source to the shock front along various lines of sight: to the top of the bubble, z\; along 
z = 0, f(0); and to the bottom of the bubble, — z^. Again, i? s , is shown with a dashed line. Together, (a) 
and (b) reveal the important feature y ~ r(0) ~ R s . (c) Thermal energy E t ^ within the bubble. Dashed line 
shows the time-evolution of £th in the spherical solution for a uniform atmosphere, (d) Pressure P within 
the bubble. Dashed line shows the solution for a spherical bubble, (e) Velocity of the top of the bubble, 
dzi/dt. The dotted line shows the internal sound speed c S; i for the parameters no — 1 cm~ 3 , Lq = 3 x 10 37 
ergs s -1 , and H = 25 pc. The Kompaneets approximation is expected to break down when dzi/dt > c s ,i 
(t > 4.9). (f) dzi/dt plotted versus z\. Dotted line as in (e). 



Table 1. UNITS OF PHYSICAL QUANTITIES 



Physical Quantity 



Unit 



Length 
Time 



(p H 5 /L y/ 3 

(Lo/poff 2 ) 1/3 



H 



Velocity 

Acceleration 

Mass 



Energy 
Pressure 



